function P_SANA=Pclogit_SANA(Y_SANA,X_SANA,Z_SANA,b_SANA)

choice=1:13;
ns=numel(choice);

Y=sortrows(Y_SANA,[1,2]);
Y=Y(:,3:end);

X=sortrows(X_SANA,[1,2]);
X=[ones(size(X,1),1),X(:,3:end)];

Z_SANA=sortrows(Z_SANA,[1,2,3]);
Z=zeros(size(Y,1),size(Z_SANA,2),ns);
for i=4:size(Z_SANA,2)
    Z(:,i,:)=reshape(Z_SANA(:,i),ns,[])';
end
Z(:,1:3,:)=[];

baseAlt=1;

P_SANA = pclogit(b_SANA,Y,X,Z,baseAlt);














